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Abstract: The unified description of fragmentation function evolution from large to small x which was 
developed for the vacuum in previous publications is now generalized to the medium, and is studied 
for the case in which the complete contribution from the largest class of soft gluon logarithms, the 
double logarithms, are accounted for and with the fixed order part calculated to leading order. In this 
approach it proves possible to choose the remaining degrees of freedom related to the medium such 
that the distribution of produced hadrons is suppressed at large momenta while the production of soft 
radiation-induced charged hadrons at small momenta is enhanced, in agreement with experiment. Just 
as for the vacuum, our approach does not require further assumptions concerning fragmentation and is 
more complete than previous computations of evolution in the medium. 
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1 Introduction 



The current description of the inclusive production of single hadrons in e + e~, ep and hadron collision 
experiments is provided by the partem model of perturbative QCD (pQCD) involving vacuum fragmen- 
tation functions (FFs) D%(x, Q 2 ) (D is a vector containing all quark FFs D q , all antiquark FFs Dq and 
the gluon FF D g ), each of which corresponds at lowest order to the probability for a parton a produced 
in the vacuum at short distance 1/Q to form a jet that includes the hadron h carrying a fraction x of the 
longitudinal momentum of a. Different theoretical approximations to the Q 2 evolution have been derived 
depending on the kinematical region of x: fixed order (FO) calculations at intermediate and large x and 
resummation to all orders of soft gluon logarithms (SGLs) at small x. In Ref. [1], an approach which 
unifies the double logarithmic approximation (DLA) at small x and the leading order (LO) Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) [2] evolution of FFs at large x was introduced. This approach 
resums both SGLs and FO logarithms in a consistent way and is defined to any order, and thus allows for 
a determination of quark and gluon FFs over a wider range of the data than previously achieved. It was 
shown to describe well both the general features of the "hump-backed plateau" [3] at small x, without 
affecting the quality of the excellent FO description of the large-x region. 

In this paper we are concerned with the more general case of medium-modified FFs in heavy-ion collision 
experiments at all possible values of x. In contrast to fragmentation in a vacuum, the fragmenting partons 
in these reactions are believed to first propagate through a quark-gluon plasma (QGP), which affects 
the FF evolution. Indeed, recent experiments at the BNL Relativistic Heavy Ion Collider (RHIC) have 
revealed the phenomenon of a strong suppression of hadrons at high transverse momentum [4, 5], in 
contrast to the expectations of the parton model and the observations from proton-(anti)proton collisions, 
which supports the picture that hard partons going through dense matter suffer a significant energy loss 
prior to hadronization in the vacuum (for a recent review see Ref. [6]). Thus the study of medium- 
modified FFs may lead to a better understanding of the physical origin of the energy loss effect by 
serving as a QGP thermometer in the nuclear matter [7]. The study of parton energy loss and medium- 
modified observables would ideally require the re-construction of jets in heavy-ion collisions. Of course, 
the huge background makes this task highly delicate. Nevertheless, thanks in particular to important 
theoretical developments on the jet re-construction algorithms [8] in a high-multiplicity environment, 
future analyses at the LHC by ALICE [9] and CMS [10] look very promising. 

Predictions for multiparticle production in such reactions at small x can be carried out by using a QCD- 
inspired model for the effect of the medium on the single particle inclusive distribution inside high 
energy jets [3], the so-called "distorted hump-backed plateau". For example, such a model was applied 
in Ref. [1 1] in the framework of the modified leading logarithmic approximation (MLLA) . In this model, 
soft emissions, which manifest themselves in the DGLAP splitting functions [2] at small x as SGLs in 
general and as 1/x infra-red singularities at LO, are enhanced at LO by multiplication of a nuclear 
medium factor N s > 1. This enhancement can be interpreted as the phenomenological description of 
the experimentally observed softening of jet spectra in nucleus-nucleus collisions. From the theoretical 
point of view, the enhancement can be seen as arising from some effective lagrangian which would be 
suitable for processes in a dense nuclear environment. 

If x is too small, the 1/x rise of the LO splitting functions is not physically correct since the FO approach 
becomes invalid. In physical applications, the effect of the enhancement on the complete resummed 
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contribution of double logarithms (DLs) to all orders, of which the l/x term is only the LO contribution, 
must be calculated by applying this enhancement to an evolution equation suitable for small x, such as 
the double logarithm approximation (DLA) that follows from angular ordering (AO) [3], or the MLLA 
master equation as was done in Ref. [11]. 

In this paper, we extend the approach of Ref. [1] to the medium in order to have a formalism which is 
suitable from large- to small- x values. This approach in the vacuum is more complete than the MLLA 
since it also incorporates the FO part of the splitting functions, which are important for describing the 
evolution at large x, and thus is more suitable for global fits. In particular, this approach allows for the 
use of small- x data to add further constraints on FFs at large x to those already provided by large-x data. 
To apply this approach to the medium involves simply repeating the same steps but in the context of the 
medium: we resum the medium-modified DLs, being the largest SGLs, in the standard DGLAP evolution 
by using the medium-modified DLA, in order to improve the small-x description of medium-modified 
evolution, while the effect of the medium on the FO part is imposed to ensure the correct description at 
large x. Just as for the vacuum, this approach can be extended to higher orders in the FO part and to 
higher classes of SGLs. Note that we do not impose momentum conservation, since partons can lose 
energy into the medium, which acts on the fragmenting partons as an external colour field. 

The rest of the paper is organized as follows. In section |2j we discuss the modifications to the FO 
splitting functions due to medium effects. Using these results, we resum the DLs in section [3] using the 
DLA equation to obtain a DGLAP evolution valid from large to small x. Finally, in section [4] we present 
our conclusions. The medium-modified form of the MLLA is derived in appendix IA1 



2 Medium-modified fixed order DGLAP evolution equations 

In both the vacuum and the medium, the DGLAP evolution takes the form 

d 



:D(x,Q 2 )= [ ld -^P{ z ,a s (Q 2 ))D( X -,Q 2 ), 

J x z \ Z / 



(1) 



where P is the matrix of splitting functions. At large x, it is accurately calculated from the FO approach 
in perturbation theory, which yields a truncation of the series 



P(x,a s ) = Y,^p( n - 1) (x) (2) 

71=1 

at some order in a s = a s j (2ir), which has the LO Q 2 evolution a s {Q 2 ) = l/(Po ln(Q 2 /Aq CD )), where 
Po = (11/6)Ca — (2/3)Tjinf is the first coefficient of the beta function and Aqcd is the asymptotic 
scale parameter of QCD. In Mellin space, defined as the integral transformation 

F(u) = [ dxx w F(x) (3) 



o 



for any function F(x), the convolution in Eq. (Q} becomes the simple product 

D{u, Q 2 ) = P(u, a s (Q 2 ))D(u, Q 2 ). (4) 



dlnQ 2 

The Mellin transform is invertible, with the inversion given by 



F(x) = ^~ [ dux-^Fiu), (5) 
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where the integration contour C in the complex oj plane is parallel to the imaginary axis and to the right 
of all singularities of the integrand F(co). By using charge conjugation and SU(ray) flavour symmetry 
to decompose the DGLAP equation into 3 simpler types of equations, Eq. (0]) can be solved analytically 
up to the desired order in a s , making the calculation of evolved FFs via the inverse Mellin transform 
numerically more efficient than the direct integration of Eq. (Q]). These 3 types of equations are the 
DGLAP equations for the non-singlet and valence quark FFs in which P is a 1 x 1 matrix, and the DGLAP 
equation for D = (D%, D g ), where P(x, a s ) is the 2x2 matrix which we write as 

(6) 

We begin with a study of the medium modifications to the latter DGLAP equation. After multiple 
rescattering in a dense nuclear medium, a relativistic partem loses a significant fraction of its energy 
scale. In the soft multiple momentum transfer model of the medium [12], the singular terms oc 1/x (as 
x — > 0) and, by the LO symmetry of the parton splitting process for all x ^ 1 under x — > 1 — x, also the 
singular terms oc 1/(1 — x) (as x — ► 1) are simply multiplied by a nuclear factor N s [11] such that soft 
emission is enhanced. Applying these modifications to the LO splitting functions in the vacuum [2,3], 
the LO splitting functions in the medium using the notation of Eq. Q become 




P^( x )=4C F (^ + fr g (x)y pj®(x) = 2C F 



1 — x 



+ /ee(x) , (7) 



p£\x) = 2C A (^ + 



x 



1 — x 



+ P^(x)=n f T R f g x(x), (8) 



where Cf and Ca are respectively the casimirs of the fundamental and adjoint representations of the 
QCD color gauge group SU(3) C , Tr = 1/2 and n/ is the number of active (anti)quark flavours. The 
plus distribution applied to a function F(x), written is defined by 



L 



i ,i 

dx[F{x)) + g(x) = / dxF(x) [g(x) - g(l)} (9) 
J 



for any function g{x). The functions f a b(x) are regular in the limit x — > 0. In particular, the Borghini- 
Wiedemann choice [11] is recovered by choosing the f a b{x) to be equal to their vacuum forms [f a b( x )]^ =\ 
given for x ^ 1 by 

x x 1 

ihg(x)] Ns=1 = [fcx(x)] Ns=1 = -- - -, (10a) 

[/<^)k=i = < l ~ x ) ~ 2 > [fMx)] Ns =i = + C 1 " x ?- ( 10b ) 

(The delta functions at x = 1 along the diagonal of the splitting functions can be determined by momen- 
tum conservation.) 

However, the correct model for the medium is not known at present, therefore neither the value of N s nor 
the fab(x) are known. We note that a more complete form for the medium-modified splitting functions 
P(x, a s ) may be derived from some effective lagrangian L me d oc (F^med) 2 m tne medium. 
The DGLAP evolution of the non-singlet and valence quark FFs is accurately described for large- to 
small-x values when the splitting function is calculated in the FO approach, since it is free of SGLs to 
all orders. Furthermore, the splitting function behaves like ~ 1/(1 — x) to all orders. Therefore, the 
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effect of the medium on the large-x LO evolution of these FFs is incorporated simply by multiplying the 
1/(1 — x) singularities in the non-singlet and valence quark splitting functions by N s , and replacing the 
rest of these functions by an unknown function which is non singular as x — ► 0, 1 (apart from possible 
delta functions at x = 1). 



In Mellin space, 



P§ g \u)=AC F \-^ + fr g {u)\, (11) 

P®(u) = 2C F (-N s Sx(u) + /se(w)) , (12) 
P$(u) = 2C A - NsS^u) + f gg (u)j , (13) 
P$(u)=n f T R f g x(u)), (14) 

where the harmonic sum Si(lo) = Yln=i V n f° r integer u>. The Mellin transform of the small- x singu- 
larity 1/x is proportional to the small-cj singularity 1/lo (note 5i(0) = 0), while the Mellin transform 
of the large-x singularity 1/(1 — x) grows at large \u)\ as lnu; (because S\{uj) ~ In a; here). Because 
P^J (u) and Pg%(oj) remain finite as \u\ approaches infinity, while 

P^H « -2C a N s Into, lu » 1, (15) 

where 

C s = C F , C g = C A , (16) 

one finds, after using Eqs. dU) and ©, the following behavior of the quark singlet and gluon FFs at large 
Q: 

D a (x,Q 2 ) ~ (1 - x) -l+4C^ s A 5 (a s (Q 2 )) ; Ae(fls) = 1 lnffls . 

The non-singlet and valence quark FFs will exhibit the same behaviour with C a = Cp- Therefore, FFs 
at large x decrease when N s increases (provided that the exponent in Eq. (fTTT ) is positive, which is the 
case at large Q). Consequently, the production of hard hadrons gets restricted at large x or, equivalently, 
in the large-p^ region in particle production in heavy-ion collisions. 

The small-x behavior of FFs is controlled by the behavior of the splitting functions in Mellin space near 
uj = 0. In the next section we will show that, in this region, the approximate behavior of D a (x, Q 2 ) as a 
function of t = In (1/x) is given by 



xD a (x, Q 2 ) ~ AAexp 
where the average multiplicity increase is given by [13] 



1 

2^' 



(18) 



Hence the production of soft hadrons gets enhanced in the medium at small x. 
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3 Evolution equations from large to small x 



In the medium, the DLA-improved evolution equation for D = D g ), which was introduced for the 
vacuum in Ref. [1] and which is valid from large to small x, becomes 

d _, _o, [ l dz'2C A N ; 



+ / 1 ^P( 2 ',o.(Q 2 ))Zj(4,Q 2 ). 12(1, 



where the matrix yl is given by 
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Equation d20l is obtained by multiplying by iV s the 1/x factor in the first term on the right hand side. 
At LO, the 1/(1 — x) singularities in the function P must be similarly enhanced, as dictated by the soft 
multiple momentum transfer model of the medium. Setting N s = 1 in Eq. (1201 ). we recover the DLA- 
improved vacuum evolution equation of Ref. [1]. Since P is free of DLs to all orders, Eq. (1201 ) can be 
used to determine the complete DL contribution to P. We note that P contains other SGLs beyond LO, 
which are less significant than DLs but significant nonetheless, and the effect of the medium on these 
SGLs is not known. 

By using the approach of Ref. [1], we will use Eq. d2"0l to obtain an evolution which is valid to DL 
accuracy at small x and to LO accuracy at large x. This goal can be achieved by repeating precisely the 
steps of Ref. [1], which we briefly outline here: We use Eq. (0) to put Eq. (1201 in Mellin space, 

where we have set D = D(u, Q 2 ), a s = a s (Q 2 ) and P = P(u>, a s ) for brevity. Substituting Eq. © into 
Eq. d2TT l and setting 

p = _pDL + pBt ( 22) 

where P DL represents the complete DL contribution to P and the remainder p^ everything else, and 
then taking the leading term (see Refs. [1, 14] for a more detailed discussion), one gets the equation 

2 (P DL ) 2 + uP DL - 2C A N s a s A = 0, (23) 

whose only solution consistent with perturbation theory is [1] 

P dl (uj, a B ) = -j (-w + V" 2 + l6C A N s a s ) , (24) 
which is the Mellin transform of 

xln± V x ) 
where J\ is the Bessel function of the first kind. 

Our unified approach to medium-modified evolution in the case of DL accuracy at small x and LO 
accuracy at large x can now be formulated, as follows: The evolution is performed using Eq. (OQ) with P 
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approximated as in Eq. G3, where P DL is given by Eq. d25l > and p^ is set equal to the LO splitting 
function after its DLs have been subtracted to prevent double counting, viz. 

P(u,a s ) = P DL (uj,a s ) + a s P^ \Lo), (26) 

where a s P (°\0) is the LO term in the expansion of P^(a;,a s ) in a s . Higher order terms in 
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N=1.8 




Inl/z 



Inl/z 



Figure 1: Quark singlet (left) and gluon (right) FFs evolved in the vacuum (N s = 1) and in the medium 
(N s = 1.8) from a factorization scale Q = 10 to 100 GeV and with the choice f a b = N s [f a b] Ns=1 . 
At the initial scale, the choice D g = (Ca/(2Cf))D-£ = iVexpf— cm 2 z]x a (l — x) 13 is made where, 
motivated by the results of Ref. [1], we have taken N = 1.6, c = 0.35, a = —2.63 and (3 = 3. 




3 lni/z 



Figure 2: As in Fig.Q] but with f ab = [f a b] Ns =i- 
P^(u;,a s ) may also be included, provided the SGLs contained within them are subtracted [1] (or 
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Figure 3: As in Fig.Q] but with = 0. 

resummed, which is not possible at present). In other words, the splitting functions in Mellin space are 
given by 

p S )(w) = w~ ( _w + ^"° 2 + 16CaNsC1s ) + 4C ^/s 9 (^), 

^IeM = -2C F N 8 S x {u) + 2C F fa{u), 

P i?H = ^ + \Au 2 + 16C A iV s a s ) - 2C A iV s £i(a;) + 2C A f gg (u), 
P$(u)=nfT R f g v(u). (27) 

The results in Eq. (1271) are the main results of this paper. 

In order to understand better the / a &, we plot the evolved FFs in the vacuum and in the medium for 
different choices of them to see which give the most physically expected results, namely an enhancement 
of hadron production at small x due to induced soft gluon radiation by the medium, and a suppression 
of this production at large x due to parton energy loss in the medium. We choose a fragmentation scale 
of 100 GeV which is suitable for the LHC. To account for medium effects, we adopt the value N s = 1.8 
from Ref. [11]. We choose the vacuum and medium FFs to be equal at 10 GeV, which is consistent 
with a finite-sized medium. In Fig. [T] we choose / a & = N s [f a t,] N =1 for the evolution in the medium 
(we absorb the delta functions at x = 1 into the definition of the f a ^ here and in what follows), i.e. the 
magnitudes of the f a b are chosen to be greater than their vacuum values. In Fig. |2j we choose the vacuum 
values exactly, f a b = [fabljy =i» while in Fig. [3] we set f a b = 0. Therefore, physically expected results 
can be obtained when the / j are chosen to not increase with N s . Note that for all choices of the /„& 
there is always a suppression at large x. 

Finally, we show the non-singlet (or valence) quark evolution in Fig. 01 Such an FF gives a measure of 
the excess of e.g. charged over neutral kaon production [15], or of positively charged over negatively 
charged hadrons. In this case there are no soft gluon logarithms at small x, and therefore there is no 
Gaussian behaviour at small x, but rather an x a behaviour, and the FO approach for the evolution is 
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Figure 4: As in Figs. [U(top left), |2](top right) and [3] (bottom), but for the non-singlet quark FF with 

N = 500, c = 0, a = 3 and (5 = 3. 

expected to be reliable from the largest- to smallest-x values. Just as for the quark singlet and gluon FFs, 
we find that the medium induces a suppression at large x and an enhancement at small x provided that 
the fat, do not increase with N s . 

3.1 Small-a; limit 

We now derive the small-x behaviour given by Eqs. ( fT8l) and ( fT9l) . The medium modifications to the 
MLLA will be derived in appendix |A] as an extension to the results below. In the small-x region we may 
neglect the second term in Eq. (l26l) . i.e. we may take 

P(uj,a s ) = P DL (uj,a s ). (28) 
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Then the solution to the gluon component of Eq. (|4TT) reduces to the simple form 



D g (u>, Q 2 ) = E gg (io, a s (Q 2 ),a s {Q 2 ))D g (io, Q 2 C 



with 



where 



E gg {u,a s (Q 2 ),a s {Ql)) = exp 



dlnQ ,2 P gg (uj,a s (Q' 2 )) 



Qo 



P gg {u, a s ) = -(-oj + y/w 2 + 16C A N S 



(29) 
(30) 

(31) 



At small x, one can predict the shape and normalization of the medium-modified FFs using the same 
procedure as for the vacuum FFs in Ref. [16], where the expression for [D a (x, Q 2 )] N _ 1 was derived 
by writing the solution in the form of a distorted Gaussian, 



xD g (x, Q 2 



Af 



cxp 



-k --sS--(2 + k)5 2 + -sS 3 + ^-k5 4 
8 2 4 V ; 6 24 



(32) 



(Tv 2vr 

In this expression, 5 = (£ — l)/cr, M is the asymptotic average multiplicity inside a gluon jet (Eq. (fT9l), 
I, a, s and k are respectively the position of the peak, the dispersion, the skewness and the kurtosis of the 
distribution. They are given by 



£ 
a 



where K n is the nth moment of D Q 



K n (Q 2 



a 3 ' 
a 4 ' 



(33) 



The reason why Eq. ( |32[ > is valid around the position of the peak (i.e. for S = 0(1)) is because the 
exponent is believed to be an expansion in a s , and this is the case in the MLLA. We now determine the 
distorted Gaussian parameters as a function of a s and also iV,. The extension to the MLLA is given in 
appendix lAl According to Eq. (|29l , the K n (Q 2 ) evolve as 

K n {Q 2 ) = AK n (a s (Q 2 ) ias (Q 2 )) + K n (Q 2 ), (34) 

where AK n (a s (Q 2 ), a s (Q 2 ) )) is the nth moment of E gg (u), a s (Q 2 ), a s (Qf ) )), given by 

AK n (a s (Q 2 ),a s (Q 2 )) = dlnQ' 2 ^-—J P gg {u, a s {Q 12 )). 

We now make the approximation that any constants and Qo -dependent terms in Eq. (l34l can be neglected, 
which is valid for sufficiently large Q. Making use of Eq. (|3"TI ). one finds 



K 2 
K 3 



SI 



-3/2 



3a. 



-5/2 



10/5 (16C A iV s )3/ 2, 



(35) 
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Then from Eqs. d33l and d35l >, one finds 



I 



4/3 



(36a) 



a 



Vm(wc A N s y/* 
o, 




(36b) 



•s 



(36c) 



k 




5 yv \ 16C A N S ' 



(36d) 



From these results, we find that the distribution around the position of the peak (i.e. 5 = O(l)) can be 
approximated by a distorted Gaussian, namely Eq. (fT8T ). The width (Eq. (I36bl >) and the kurtosis (Eq. 
d36d| |) of this distorted Gaussian are suppressed in the limit N s 3> 1 by 1/Nf, where a is a positive 
constant. 

Setting oj = in Eq. (|3TT >. one recovers the rate of the mean average multiplicity increase in a gluon jet 
with respect to In Q 2 , i.e. the DLA medium-modified anomalous dimension of the multiplicity [13]: 



in which case Eqs. d29l ) and (l30l) give Eq. ( fT9l ) at large Q. 

4 Conclusions 

Our approach provides a general framework for the incorporation of medium effects into the DGLAP 
evolution of FFs from large to small x, which will be important for the description of single hadron 
production in heavy-ion collisions from high to low pr- In particular, our approach allows the use of 
measurements of such observables at low pt, e.g. from heavy-ion collisions at the LHC, to provide 
additional constraints on the medium-modified FFs to those constraints provided by large-pj- data, and 
therefore could improve the ruling out of various models of the medium by improving the constraints on 
the unknown degrees of freedom. 



The MLLA can be regarded as a simplification to the approach of Ref. [1] and of this paper, in which the 
main qualitative features of hadron production are preserved. We now use our approach to determine the 
form of the MLLA with medium effects taken into account. The MLLA is obtained by setting u> = in 
the second term in Eq. (|26"T ). i.e. by taking 



P gg (0,a s ) = ^/CA~N ^ 



(37) 



A MLLA limit 



P(u,a s ) 



P DL ( W ,a s )+a s P^(°)(0) 



(38) 



which for small uj is correct to terms of 0(ui). In the case of the medium, 




n f T R f g x(0) 2C A f gg (0) 




(39) 
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which corresponds to the hard SL corrections in the MLLA formalism [3]. Setting 



the system in Eq. (121]) with Eq. (l39l > reduces to a self-contained equation for D g and a coupled equation 
for constrained by the solution of D g , which can be written in the form 

W + 2 dl^) = ^ F N a a s D 9+ Uc FhM + 2 ^fM^ + 2^1^,(40) 



(41) 



In Ref. [11] the choice N s = 1 is made in P BC (°)(0), which is known: 
and in a mec i, which is also known: 



However, we will keep N s / 1 in what follows in order to maintain generality. Using Eqs. 
Eq. (PiTI ) after some straightforward algebraic operations becomes Q 



dP 1 

2f3(a s )-^ - -a med {uj + 2P gg )a s 



da s 2 

which provides the following solution correct to terms of single logarithmic order: 



P gg (uj, a s ) = - (-U) + \/uj 2 + lQC A N s a t 



+ GU 



AC A N s a s 



2 + IQCaN s 



UJ 



1 / UJ 
Turned 1 H ; = 



(42) 



and(|30b. 



(43) 



(44) 



Setting uj = in Eq. (1441) . one recovers the rate of the mean average multiplicity increase in a gluon jet 
with respect to In Q 2 , i.e. the MLLA medium-modified anomalous dimension of the multiplicity [13]: 



Pgg(0, Os) = V CAN s a s + (f3 - a me d) • 



(45) 



The purpose of (3q is to account for the running of the coupling constant a s . Finally, the solution for D% 
at NLO in Mellin space can be derived from Eq. d40l ), once D g has been computed, and reads 



2C F f 1 



1 

^T-amed + /Sg(0) + ^-/es(0) 



uj + 2 



d 



dlnQ 2 



(46) 



Since [uj + 2 



'dlnQ* 



~ 0(\/N s a s ), the second term on the right-hand side of Eq. (l46l ) is of order 



0{yf a s /N s ). For N s = 1, one recovers the result in the vacuum [17] with 



1 C F 
^r-amed + /s fl (0) + ^-/ss(0) 



iV s =i 



4CU 



'A term +a mec |/3(a s ) on the right hand side has been neglected since it is not required in what follows. 
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Making use of Eq. (I44l ). one finds 



-1 -1/2 
7^ 0, s a med a s 

K l = T77" + 



4/3 2A)V16C^JV S ' 

3 /2 „_i 



(Is 



K; 



-3/2 



2/3 (16C A ^)3/2' 

3a s 5 ^ 2 3a~ 2 
10/? (16C A iV s )3/ 2 + (WC A N s f 



K4 = ~ — 7-77; » T vWo + ,n^\ M2 - ( 4? ) 



Then from Eqs. d33T > and (|47T ). one finds 



N = No exp 



2 C A iV 1 

+ 777- (flmed - Po)lna s 



" - V!^(i6C A jv,)V4 ^ " 2*V iscM 1 ' (48c) 



3/4 



(48a) 



/3 V a s 4/? 

-1 -1/2 
^- + ame * as (48b) 

-3/4 / o 

a. ' f. 3 



-^"-lieckl ' (48d) 



which coincide with the results in Ref. [16] in the limit N s = 1. 

Assuming that a me( j rises slower than ivf //4 , the skewness (Eq. (I48dl >). like the width (Eq. (I48cl l) and the 
kurtosis (Eq. (I48el >). is suppressed in the limit N s S> 1. As found in subsection 13. II the distribution close 
to the position of the peak can be approximated by a Gaussian shape. Finally, in Ref. [11], the choice 
[flmed]7v s =i = a * s ma de, in which case the position of the peak approaches the asymptotic DLA value 
%fofar large N s . 
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